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The pulse displacements produced by timing errors in a chain of regen- 
erative repeaters (using tuned-circuit timing filters) are represented by a 
linear transformation of the pulse displacements at the output of the first 
repeater. To facilitate the discussion of the general problem, the simpler 
case of periodic pulse trains is considered first. For this case it is shown 
that while the mean value tends to infinity, the central moments of the pulse 
displacements remain bounded as the number of repeaters approaches 
infinity. Further results are obtained which show that all the moments of 
the spacing jitter remain bounded for an indefinitely long string of re- 
peaters. Finally, the misalignment in the jitter at any given repeater is 
represented by a simple expression which shoivs that the essential com- 
ponent in the misalignment is flat delay. 

The general problem of random pulse trains, infinite in length, is dis- 
cussed in Part II in this issue. The residts obtained for the general case 
are quite different from those obtained for the periodic case. The variance 
is unbounded in this case except for pulse trains with certain special re- 
strictions. The computational aspects for the evaluation of jitter accumula- 
tion will be discussed in a subsequent paper. 

I. INTRODUCTION* 

In regenerative digital transmission systems, one of the important 
problems is that of maintaining the proper distance between the signal 
pulses. The problem becomes much more serious when the system con- 
tains a rather long chain of regenerative repeaters. Several aspects of a 
theoretical nature in connection with this problem have been discussed 
by Sunde, 1 Bennett, 2 Rowe 3 and Rice. 4 

We study here the pulse displacements produced by tuning errors in a 
chain of repeaters using tuned-circuit timing filters. For simplicity, we 
shall consider the system free of noise, distortion, etc. 
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An idealized version of the physical system is a chain of repeaters 
with the input supposed to be a train of unit impulses. Each repeater 
is a device containing a resonant circuit which is excited by the incoming 
train of pulses. The response of the resonant circuit to the incoming 
signal will ideally consist of a sum of sinusoids and will pass upwards 
through zero at an instant determined by the resonant frequency of the 
circuit. This instant will coincide with the instant of occurrence of the 
pulse, if it occurs at all, when the resonant frequency is identical with 
the pulse repetition frequency. The repeater does its "repeating" by 
sending out a unit impulse, at the instant the response of the tuned- 
circuit passes upwards through zero, provided the input signal has a 
pulse at or near the same instant. If there is no pulse in the input, no 
pulse is sent out. 

Due to tuning error, the tuned circuit in a practical repeater would 
resonate at a frequency somewhat different from the pulse repetition 
frequency. Further, the impulse response of the circuit is more truly a 
damped sinusoid. These considerations show that the positions of the 
pulses sent out by a practical repeater are somewhat displaced from the 
true positions of the pulses in the original pulse train. 

Actually the system consists of a chain of repeaters. We are thus led 
to a consideration of the statistical properties of the pulse displacements 
produced in a random pulse train by the combined effect of mistuning 
in each successive repeater. Of particular concern is the behavior of 
the pulse displacements as the number of repeaters gets larger and larger. 
It is to this question that we attend. 

We begin our discussion by a mathematical statement of the problem. 
We show that the pulse displacements at the output of a chain of 
repeaters may be represented by a linear transformation, in a Banach 
space, of the pulse displacements at the output of the first repeater. 

The linear operator (or, the linear transformation) becomes un- 
bounded, in the limit, as the number of repeaters gets indefinitely large. 
From this follows the result that the average value* of the pulse dis- 
placements increases indefinitely as the number of repeaters approaches 

infinity. 

The behavior of the variance, as well as the other central moments, 
of the pulse displacements is investigated by considering a suitable 
projection, when it exists, in the Banach space. When the domain of 
the above linear transformation is a linear manifold obtained by the 
desired projection, we find that the linear operator is bounded. Conse- 

* All averages are taken over the values of the pulse displacements. No averages 
over the mistunings should be compared with the results obtained here. 
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quently, all the central moments of the pulse displacements are shown 
to remain boimded as the number of repeaters approaches infinity. 
When the above-mentioned projection does not exist, the central mo- 
ments are shown to be unbounded. 

Practical situations call for a determination of the bounds on the 
central moments when the number of repeaters is finite. In such cases, 
the input pulse trains may be assumed to be periodic pulse trains with 
the period much larger than the time constants of the timing filters. 
The problem reduces to a linear transformation in a finite dimensional 
vector space. The central moments are bounded and they can be pre- 
cisely evaluated. A simple procedure to determine these bounds is 
developed. 

The same analysis can be directly applied to an investigation of the 
so-called "spacing jitter," or variations in the spacings between virtual 
pulse positions. Similar results are obtained for both a finite and an 
infinite number of repeaters in the chain. 

We shall also have occasion to remark upon the "misalignment noise" 
which is the jitter introduced, by the nth repeater, in an already jittered 
pulse train coming into the same repeater. 

Finally, in a subsequent paper we shall discuss the computational as- 
pects for the evaluation of jitter accumulation in a long string of re- 
peaters. 

Organization of the paper is as follows. We start with the statement 
of the problem in completely general terms and express it as a linear 
transformation. Next, to facilitate the discussion of the general prob- 
lem, we consider the simpler case of a periodic pulse train. In Part II of 
the paper,* we consider the general case of a completely random pulse 
train. 

II. STATEMENT OF PROBLEM 

The input to the chain of repeaters is supposed to be a train of unit 
impulses which occur, if they occur at all, at the instants 



{■•-,-27,-7,0, 



t. it. 



The occurrence or nonoccurrence of a pulse at time t = — m is deter- 
mined by the value of the random variable Ct„ . If G„ = 1 , which happens 
with a given probability, a pulse is present. If G„ = 0, no pulse is pres- 
ent. 

The resonant circuit in the repeater is excited by the incoming train 

* Part II of the paper appears in this issue, p. 1781. 
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of pulses. The response of the circuit to a unit impulse at time t = is 
assumed to be 

e~' 1 sin cod, I > (1) 

where w r is almost 2w but, due to tuning error, misses its desired value 

by 

« = 2ir - u„r = 2ir(./V - h)/j r . (2) 

Here /, = 1/t is the pulse repetition frequency. The decrement a is 
related to the Q of the circuit by <jt = r/Q. 

The response of the resonant circuit to the incoming pulse train will 
consist of a sum of terms of the form ( 1 ) and will pass upwards through 
zero at an instant near t = (—tit), say at t = (—/it + L). The re- 
peater sends out a unit impulse at the instant (—tit + t„) if the input 
signal has a pulse near (— nr). If there is no pulse in the input, no pulse 
is sent out. The response of the resonant circuit still goes through zero, 
and we can say that there is a "virtual" pulse displacement of amount 
t„ seconds (or of 2irt„/T radians). 

For a chain of repeaters, we assume that all of the resonant circuits 
have the same Q but that their mistunings « , e 2 , • • • are distributed 
independently and at random. Let ft' be the displacement of the fcth 
pulse (originally entering the first repeater at t = —kr) as it comes out 
of the Zth repeater where / = 1, 2, • • • . The displacement &' is measured 
in radians, where 2* radians corresponds to the pulse interval t. The 
superscript I signifies the output of the /th repeater. The mistiming in 
the resonant circuit in the Ith repeater is represented by ej . When we 
assume that Q is very large and the mistunings e t are much smaller than 
aT = T /Q radians, we are led to a set of equations which relate the pulse 
displacements out of the Ith repeater to those out of the (I - l)th re- 
peater. These equations are 

00 
ft — 



22 <W (3) 

ii -0 

(I = 1,2,3, •••;* -0,1,2, ■••), 

where = exp ( — <rr) « 1 - (*■/<?) is a number slightly less than unity. 
The initial conditions are that the pulses entering the first repeater 
have zero displacement, i.e., 

g fc ° = 0, fc - 0,1,2, .... (4) 
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These equations are given by Rowe and also by Rice. 4 Here we have 
followed their terminology very closely. 

The physical problem dealing with a chain of repeaters is now re- 
placed by the mathematical problem of studying the behavior of the 
variables £ t defined by the above equations. The (2,,'s and ej's are either 
given explicitly or are random variables whose distributions are known. 

III. LINEAR TRANSFORMATIONS 



We note that the set of equations in (3) is a linear set, and we can 
express it as a linear transformation of the set of variables |&' -1 ) into 
the set {& } . We are, however, primarily interested in the behavior of 
{&-'} when / is large and when no knowledge of (& l ] is available. A 
more useful expression is obtained by rewriting (3) as 



6/- 



E G-+*0"$» 



+k 



E G«+*/3" 






(5) 



where 



& — 



2 Gn+kP"nei 



Z «n+t/S" 



(6) 



In our formulation, zero mistiming does not introduce any jitter in a 
jitter-free pulse train. We will therefore understand the chain to start 
with a repeater having non-zero mistuning. 

Equation (5) can be used to express {£*.'} as a linear transformation 
of !& J. To do this, define a matrix (infinite) 
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where 



1774 THE BELL SYSTEM TECHNICAL JOURNAL, NOVEMBER 1962 

Si = Z ««+.£"; (8) 

71 = 

and define a vector 

x t = &UU', •••]• (9) 

Then (5) becomes 

X t = TX l _ 1 + € !-X 1 ; (Xo = 0), (I = 1,2,3, ■• •)• (10) 
From (10) it follows that 

Xr = \- ! Le l -X']x 1 , (T°- /). (11) 

One can, if need be, discuss the behavior of (11) in the above form. 
However, the e's are usually of the same order of magnitude, and the 
equation is considerably simplified by assuming that the e's are iden- 
tical.* Then 

x l+1 = I S r\ Xi . (12) 



-!»• 



We are interested in the problem when I becomes indefinitely large, or, 
dropping superfluous subscripts, 

Y = lim I E T I X. (13) 



Hm [i: r\ x. 

l-oo \_y=0 J 



Here X and Y represent the pulse deviations out of the first repeater 
and out of the (I + l)th repeater, respectively. 

The original problem is now represented as a linear transformation of 
X into Y. The linear transformation, when it exists, is a function of 
another linear transformation T. The domain, as well as the range, 
of the transformation T is a Banach space, as will be shown in Part II. 
Here, we pursue the simpler case of a periodic pattern. 

Whether the variance is bounded or not is not a particularly impor- 
tant question for the periodic case. Such a question can be answered 
by a very simple argument. However, we give here instead a complete 
analysis of the periodic case. Our purpose in doing so is twofold. First, 
the analysis shows how certain basic properties of the operator T in- 
fluence the questions of boundedness of the jitter; it also gives a simple 

* We shall discuss elsewhere the difference, if any, in the results when we do 
not make this assumption. 
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computational procedure for evaluating the accumulated jitter. Second, 
the analysis serves as a simple introduction to the more complex argu- 
ment pursued in Part II. 



IV. PERIODIC PULSE TRAINS 



We assume here that the a n 's form a pattern which repeats itself 
with a period m. The pattern is otherwise arbitrary. In such cases, the 
pulse displacements are also periodic with the same period m. Then, 



(14) 



a„ +m = a n 

£ft+m = £<t 

for all values of indices n and k. 

The domain of the operator T is thus an m-dimensional space. Since 
Gin+m = &n , the range of T is also of dimension m. The problem reduces 
to the study of a linear transformation in a finite dimensional space. 
The operator T is now represented by a finite matrix A. 
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where 

8*' = (1 - /3 m )s t . 
For the periodic case, (13) becomes 



Y = lim IE 
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(15) 
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where we continue to use the same symbols X and Y to represent the 
finite dimensional vectors. 

In order to investigate the limit of (17), one must determine the 
behavior of the infinite series and its convergence properties. Moreover, 
if the limit does not exist, the question to be answered is whether or not 
the variance of Y has any limit. Other central moments may also be of 
interest. 

In what follows, we show that the limit of (17) does not exist. This 
implies that the mean of Y is infinite. However, we shall show that the 
central moments always exist for any arbitrary m. We assume through- 
out this paper that averages over the sample values* are statistically 
identical to the averages over the ensembles. 

A discussion of the properties of the linear transformation defined by 
(17) involves the study of a function of the matrix A. In order to discuss 
such a function, one must have a knowledge of the spectrum of the 
matrix. We study the spectrum of A in the next section. 

V. SPECTRUM OF A 

In this section, we prove the following theorem. 
Theorem: The spectrum of A consists of two parts: 

1 . The maximum eigenvalue is located at X = 1, and it is simple; 

2. All other eigenvalues are such that their modulus is less than unity, 

i.e., | X t - 1 < 1. 

Proof: Observe that A is a stochastic matrix since the sum of each 
row is equal to one and all the elements of the matrix are nonnegative. 
Thus, X = 1 is indeed an eigenvalue with eigenvector {1, 1, • - • , 1}. It 
also follows that the entire spectrum of A is contained in the unit disk 
| X | ^ 1. This can be observed in a simple manner by considering powers 
of matrix A and noting that the trace of A" does not exceed m, the 
order of the matrix A. If there were any eigenvalue for | X | > 1, one 
could find a large enough n such that the trace of A" would exceed m. 
(We do not worry about cancellation because we can always choose 
the proper n to prevent this.) Hence, there are no eigenvalues outside 
the unit disk. 

Next, we wish to show that there are no other eigenvalues (X y* 1) 
with modulus equal to one. We obtain a matrix equivalent to A by 
means of elementary transformations of interchanging rows as well as 
the corresponding columns. The eigenvalues of the matrix are invariant 

* The values of the pulse displacements are referred to as the sample values, 
and the ensemble is the set of admissible sequences of pulse displacements. For 
justification of the above assumption in the general case, see Bennett, op. cit. 
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under such operations. We obtain a matrix of the form 



B = 



r A' 



C D 



(18) 



where A' is a square matrix all of whose elements are positive and D 
is a square null matrix. Only the eigenvalues of A' need be considered. 
To A' we apply Perron's theorem which, for a stochastic matrix with 
all elements positive, states that: the extremum eigenvalue is located 
at X = 1 ; it is simple; and its modulus exceeds the moduli of all other 
eigenvalues. Q.K.I). 

VI. MEAN, VARIANCE, ETC. 

The solution to (17) can now be expressed in terms of the basis vectors 
of A in the form 

Y = lim £ EXpVY*" 1 , O 9 ) 

where, X is the eigenvector of A corresponding to the eigenvalue 
X,, of A . The coefficients a„ are the expansion coefficients in 

X = E ^A' (M> . (20) 

We have assumed, for the present, that A is of simple structure. There 
are no significant changes in the development when such an assumption 
is not made. We shall discuss this matter a little later. 

In the previous section it has been proved that the extremum eigen- 
value, say Xi , is simple and is located at Xi = 1 . The rest of the eigen- 
values are strictly inside the unit circle. The mean value of Y is seen to 
approach infinity by considering only those terms that involve X x = 1, 

M m 

Y = ai X ll) E Xi' + £ i-^V X&, (21) 

where, A' (,) = { 1, 1, • • • , 1). 

The first term in (21) is a divergent series and }" approaches infinity 
as the number of repeaters increases indefinitely.* The behavior of the 
central moments is investigated by considering 

[Y - f ] - Z (i^ry) \X W - XWJ- (22) 



* The statement is valid, in general, provided «i ^ 0. We need only show that 
there exists at least one X such that o, j± 0. Consider a pulse train with all pulses 
present; then X = «iX (1) with a\ j* 0. 
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The value of [ Y - Y) is finite since | X„ | < 1. All powers of [Y - ?] 
are finite and [Y — Y] h will also be finite. We thus see that all the cen- 
tral moments, including the variance, are finite. 

Now we consider the case when the structure of A is not simple. The 
only difference in this case concerns the vectors corresponding to eigen- 
values other than Xi . Let us, for simplicity, consider the basis vectors 
that correspond to an eigenvalue X„ of multiplicity two. Similar develop- 
ments can be carried out when the multiplicity is greater than two. 
The normal form of A would have a Jordan block 



X 



(23) 



It is well known that there exist two linearly independent vectors 
Xi W and X^ such that 

AXl - X " Al (24) 

AX, 00 = X M X 2 <A ° + Zi <M) . 

The vector Xi W is an eigenvector of A and is transformed in the same 
manner as the vectors X w are, and it yields for Y a term of the form 

(A)*" (25) 

On the other hand, when A" operates on X 2 W it yields 

A'Xf* = X/X. 00 + *V- , X 1 ( "\ (26) 

Thus X 2 ( " ) contributes to Y a term of the form 

«/ £ A'X, W = «/ \t, X/X. w + E vxr'x^'l , 

y=0 [_V=0 »=0 J 

since | X„ | < 1. 

The terms due to the basis vectors of .4 corresponding to X„ are shown 
to be bounded, and our results on the boundedness of the central mo- 
ments remain valid regardless of the structure of A . 

VII. SPACING AND MISALIGNMENT 

Sometimes a more useful measure of jitter is the spacing jitter, which 
is defined as the deviations in the spacing between adjacent pulses or 



(27) 
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pulse positions. This is obtained by taking the difference of the ad- 
jacent pulse position deviations. To do this operationally, let us define 
an operator S which shifts the elements in the vector Y such that the 
fcth element appears as the (A: — l)th element and the first element 
appears as the with element . 



S = 



10 
10 
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10 









(28) 



The spacing jitter Y, can then be represented in terms of the timing 
jitter Y by 



Y„ = \I - S\Y. 
By using (19) and (29), we have 

l 7,1 

Y. = lim £ £ \JaJJ - S}X W . 

/-.00 y-Q ^=1 

The operator [/ — S] annihilates X and we obtain 
F. = Z (r-^r) [I ~ SIX'*- 



(29) 



(30) 



(31) 



The spacing jitter is finite for all sample values, and so the mean and 
all other moments of this jitter are finite.* 

Next, we briefly consider the misalignment which is defined as the 
difference between the timing errors at the output and at the input of 
a given repeater. The representation of the misalignment in the (I + l)th 
repeater is given, in the periodic, case, by modifying (12) to a finite 
dimensional one and obtaining 

[X t+1 - X,} = A l X lf (32) 

where, Xt represents the jitter at the output of the fcth repeater. 

Equation (32) implies that the misalignment essentially amounts to a 
flat delay as / gets larger. Indeed, there is virtually no difference in the 
misalignment for different repeaters when the values of / are reasonably 
large. 



* For periodic pulse patterns, this is intuitively obvious. 
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VIII. CONCLUSION 

The general problem of timing errors in a string of repeaters has been 
expressed in terms of certain linear operators and functions of these 
operators. The simpler case of periodic pulse patterns is then studied in 
detail. We have shown, for the periodic case, that the mean value of 
jitter accumulation in a string of repeaters increases indefinitely but 
that the central moments of the jitter remain bounded. In fact, the 
divergence of the mean value for the infinitely long string stems from 
the accumulation of the flat delay occurring in each repeater. Once this 
flat delay is eliminated, the remaining part of the jitter is bounded. 
Consequently, all the central moments are bounded. All the moments 
of the spacing jitter are bounded for identical reasons. The misalign- 
ment behavior is also explained by the dominance and the invarianee 
of the flat delay. 

The question of evaluating the jitter accumulation will be discussed 
in a subsequent paper. We will show there that the spectrum of the 
operator A can be determined fairly simply even for very large periodic- 
ity. No polynomials of high degrees need be solved to determine the 
eigenvalues. We shall also discuss the computation errors involved in 
periodic approximation versus those involved in truncation of the 
infinite pulse train. 

The general case of random pulse trains with no periodic structure 
will be examined in Part II. We shall have occasion to thoroughly 
examine the operator T. Since we shall be concerned with infinite 
dimensional space, the spectral properties of T are not so easy to de- 
termine. We shall compare the spectral properties of T with those of A 
in order to delineate the difference between the two cases. 
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